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Abstract 

We extend the earlier suggested QCD-motivated model for the (^-dependence 
of the generalized Gerasimov-Drell-Hearn (GDH) sum rule which assumes the 
smooth dependence of the structure function gj-, while the sharp dependence is 
due to the 52 contribution and is described by the elastic part of the Burkhardt- 
Cottingham sum rule. The model successfully predicts the low crossing point 
for the proton GDH integral, but is at variance with the recent very accurate 
JLAB data. We show that, at this level of accuracy, one should include the 
previously neglected radiative and power QCD corrections, as boundary values 
for the model. We stress that the GDH integral, when measured with such a 
high accuracy achieved by the recent JLAB data, is very sensitive to QCD power 
corrections. We estimate the value of these power corrections from the JLAB 
data at Q 2 ~ lGeV 2 . The inclusion of all QCD corrections leads to a good 
description of proton, neutron and deuteron data at all Q 2 . 
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The generalized (Q 2 -dependent) Gerasimov-Drell-Hearn (GDH) sum rules ^ |2] 
are just being tested experimentally for proton, neutron and deuteron [3 IH El- 
The characteristic feature of the proton data is the strong dependence on the four- 
momentum transfer Q 2 , for Q 2 < lGeV 2 , with a zero crossing for Q 2 ~ 200 — 250MeV 2 , 
which is in complete agreement with our prediction j3|Hj, published almost 10 years 
ago. Our approach is making use of the relation to the Burkhardt-Cottingham sum rule 
for the structure function g 2 , whose elastic contribution is the main source of a strong 
(5 2 -dependence, while the contribution of the other structure function, gr = 9i+ 9i is 
smooth. 

However, the recently published proton JLAB data |Hj lie below the prediction, 
displaying quite a similar shape. Such a behaviour suggests, that the reason for the 
discrepancy may be the oversimplified treatment of the QCD expressions at the bound- 
ary point Qo ~ lGeV, defined in the smooth interpolation between large Q 2 and Q 2 = 
and which serve as an input for our model. For large Q 2 we took the asymptotic value 
for the GDH integral and we neglected all the calculable corrections, as well as the 
contribution of the g 2 structure function. This was quite natural and unnecessary 10 
years ago, since no data was available at that time. 

In the present paper we fill this gap and include the radiative (logarithmic) and 
power QCD corrections. We found that the JLAB data are quite sensitive to power 
corrections and may be used for the extraction of the relevant phenomenological pa- 
rameters. We present here the numerical values of these parameters which naturally 
depend on the approximation of the QCD perturbation theory. The resulting theoret- 
ical uncertainty should be of order of the last term of the perturbative series , taken 
into account, and should not therefore be more than several percents. 

Moreover, the perturbative series should contain the renormalon ambiguity due to 
the factorial growth of the coefficients, resulting in a power rather than logarithmic 
correction with an unspecified coefficient. It is in fact this ambiguity which allows 
the interpretation of the dependence of the numerical value of the power correction, 
earlier mentioned for the case of the F3 structure function j^J, as an ambiguity in the 
separation of logarithmic and power corrections. 

We use the values of the power corrections as an input for our model at Ql ~ lGeV 2 
and we achieve a rather good description of the proton data at lower Q 2 . We also 
present the improved description of the neutron and deuteron data and the behaviour 
of the Bjorken sum rule at low Q 2 . 

The starting point of our approach is the analysis of the general tensor structure of 
Wjf, the spin-dependent part of hadronic tensor W^ v . It is a linear combination of all 
possible Lorentz-covariant tensors, which should be orthogonal to the virtual photon 
momentum q, as required by gauge invariance, and linear in the nucleon covariant 
polarization s, from a general property of the density matrix. If the nucleon has 
momentum p, we have as usual, s ■ p = and s 2 = —1. There are only two such 
tensors: the first one arises already in the Born diagram 

TT = e^s a q p (1) 

and the second tensor is just 

Tf = (s ■ q)t^p a q p . (2) 
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The scalar coefficients of these tensors are specified in a well-known way, since we have 



= qp{gx{%, Q' z )s a + g2(x, Q 2 )(s a - p a :LJL )) 

p ■ q p ■ q 



-QpiMx, Q 2 ) + g 2 {x, Q 2 ))s a - g 2 (x, Q 2 )p a ) . (3) 



p ■ q p ■ q 

Therefore, due to the factor (s ■ q), g 2 is making the difference between longitudinal 
and transverse polarizations, while gr = gi + #2 contributes equally in both cases. 
Let us consider now the Q 2 -dependent integral 

2M 2 9 /If 2 f 1 

HQ 2 ) = ^rlMQ 2 ) = f£L / 9l (x,Q 2 )dx . (4) 



Q 2 v " ' Q 2 

It is defined for all Q 2 , and gi(x, Q 2 ) is the obvious generalization for all Q 2 of the 
standard scale-invariant structure function gi(x). Note that the elastic contribution at 
x = 1 is not included in the above sum rule. Then, by changing the integration variable 
x — > Q 2 /2Mu, one recovers at Q 2 = the integral over all energies of spin-dependent 
photon-nucleon cross-section, whose value is defined by the GDH sum rule PjJ |2] 

h(0) = , (5) 

where \xa is the nucleon anomalous magnetic moment in nuclear magnetons. While 
ii(0) is always negative, its value at large Q 2 is determined by the Q 2 independent 
integral gx(x)dx, which is positive for the proton and negative for the neutron. 

The separation of the contributions of gx and g 2 leads to the decomposition of 
h{Q 2 ) as the difference between It{Q 2 ) and I 2 {Q 2 ) 



h{Q 2 ) = l T {Q 2 )-h{Q 2 ), (6) 



where 



9M 2 f 1 2M 2 f 1 

MQ 2 ) = -qT J o 9t(x, Q 2 )dx, I 2 (Q 2 ) = -q 2 -J o 92(x, Q 2 )dx . (7) 

There are solid theoretical arguments to expect a strong Q 2 -dependence of I 2 {Q 2 ). 
It is the well-known Burkhardt-Cottingham sum rule ^U]. It states that 

I 2 (Q 2 ) = ^G m (Q 2 ^ Gm(Q2) ~ Ge(Q2) 



4^ ' 1 + Q 2 /AM 2 ' v ' 

where \i is the nucleon magnetic moment, Gm{Q 2 ) and Ge{Q 2 ) denote the famil- 
iar Sachs form factors, which are dimensionless and normalized to unity at Q 2 = 0, 
Gm(0) = Ge(0) = 1. For large Q 2 , as a consequence of the Q 2 behavior of the r.h.s. 
of Eq. (jHJ) , we get 



/ g 2 (x,Q 2 )dx = 0. 
Jo 



(9) 
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In particular, from Eq. ® it follows that 



J 2 (0) = £±^, (10) 

e being the nucleon charge in elementary units. To reproduce the GDH value (see 
Eq. (jSJ)) one should have 

MO) = ^f, (ii) 

which was indeed proved by Schwinger [TJ]. The importance of the g 2 contribution can 
be seen already, since the entire /z^-term for the GDH sum rule is provided by I 2 - 

Note that It does not differ from I\ for large Q 2 due to the BC sum rule, but it is 
positive in the proton case. It is possible to obtain a smooth interpolation for T T (Q 2 ) 
between large Q 2 and Q 2 = [7] 



,» A , P 2M 2 Q 2 2M 2 

(Ql) 2 1 Q 2 



i p T (Q 2 ) = 0(Qo - Q 2 )(^r - -tt^K?) + e (Q ~ Qo)^rK (12) 



where = f g\(x)dx. The continuity of the function and of its derivative is guaran- 
teed with the choice Ql = (16M 2 / fj, A> p)T^ ~ lGeV 2 , where the integral is given by the 
world average proton data. 

This smooth interpolation seems to be very reasonable in the framework of the 
QCD sum rules method [HIE], as the low energy theorem for the quantity linear in \ia 
may, in principle, be obtained by making use of Ward identities. It is also compatible 
with resonance approaches [13 , as we observed earlier [S] that the magnetic transition 
to A(1232), being the main origin of sharp dependence in that approach, contributes 
only to 02 • 

However, such interpolation neglects the QCD perturbative and power corrections 
and on the other hand, it assumes that at the boundary point Qo the contribution of 
g 2 is already extremely small so that 

i p AQl) = mQ 2 o)- (is) 

Both types of corrections are easily taken into account, although this does not allow a 
simple analytic parametrization. 

The starting point of the upgraded model is the corrected expression for the asymp- 
totic expression for I{ (i = p,n): 



il(Q 2 ) = 2 -^\j\\M 2 )dx{i - ^ 1 )-c, << qI >> ] , (14) 

where we took into account the one-loop perturbative correction (while the inclusion 
of higher order ones will be discussed later), as well as the twist-4 contribution j!4j . 
Here q is the charge factor equal to 2/9 for proton and to 1/18 for neutron, while the 
matrix elements of the combinations of reduced twist-3 and -4 operators happen to be 
equal jHj for both proton and neutron: 

<< O p »=« O n »= 0.09 ± 0.06GeV 2 . 
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Note that the kinematical target mass corrections happen to be numerically small 
and we neglect their contribution. Anyway, they may be combined with the genuine 
twist corrections and the resulting change of the latter is within experimental and 
theoretical errors. 

As the expression for the I2 stays unchanged, the expression for I? above the 
matching point Qq should change accordingly. Let us start from the proton case. 

i p T , pe JQ 2 ) = o(Q 2 - <S)^W - ^) - ^qI^ ~ W 2 )]- (15) 

The smooth interpolation to the GDH value at Q 2 = is now more difficult and cannot 
be performed anymore, by making use of simple analytic formulae. Instead, we expand 
()15j) to the power series at the point Qq and define the expression at the low Q 2 as: 

N / 1 d n T p \ 

i p t„AQ 2 ) = o(Ql - Q 2 ) £ (Q 2 - QlT- (16) 

n=0 V U ^ > ' Q=Qo 

Here N is the number of continuous derivatives of these two expansions, which turns 
out to be a free parameter of the model, together with the matching value Qq. They 
should be chosen in such a way, that the condition for real photons 

I P T„rM = ^ (17) 

is satisfied. 

The procedure we are implementing in such a way, may be considered as a matching 
of the "twist-like" expansion in negative powers of Q 2 and the "chiral-like" expansion 
in positive powers of Q 2 , which is similar to the matching of the expansions in direct 
and inverse coupling constants. In its simplest present version we take only the value 
1(0) as an input, although the slope and other derivatives calculated within the chiral 
perturbation theory may be added in future work. 

As soon as the low Q 2 region exhibits the important contribution of the resonances 
[13], the suggested procedure may be also considered as a version of quark-hadron 
duality. It is worthy to note here, that Bloom-Gilman duality for spin-dependent case 
is strongly violated by the contribution of A(1232) resonance jT^j. As it was mentioned 
before, since this resonance does not contribute to the structure function 0, it is 
this function which may be a good candidate to study duality. 

We have studied Eq. (fTTj) numerically changing the following inputs: 

i) for different order of perturbative correction (1,2,3 loops) |lfij . 

ii) for different values of the degree of approximating polynomial N in Eq. (|16p: it 
is interesting that taking N — 1 does not allow for solution of Eq. ([17)1 . 

iii) for different values of non-perturbative corrections, which we were choosing in 
order to be close to JLAB data at their highest Q 2 ~ lGeV 2 . We observed that the 
increasing of the order of perturbative corrections lead to systematical decrease of the 
required non-perturbative one, which is similar to the case of _F 3 structure function [9 j 
and may be considered as a manifestation of the ambiguity in separating logarithmic 
and power corrections. 

iv) we varied the matching point Qq until Eq. (fTTj) is satisfied. 
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We found that Q is systematically (but not strongly) increasing with N. The 
expression for Q 2 dependent integral Y\{Q 2 ) = If non _ pert (Q 2 )Q 2 /2M 2 resulting from 3- 
loops perturbative correction with N = 3, << O p »= O.llGeV 2 and Ql = 0.97GeV 2 
is shown in Fig. 1. It is reasonably close to the JLAB data j^j. In what follows the 
thick lines correspond to our new approach and we present also the results from the 
old approach for comparison. 



. 1 




. 02 



Figure 1: Our description of Y\(Q 2 ) versus Q 2 . The thick line is the new analysis 
to be compared with the thin line, which represents our previous approach without 
corrections. 

Here we took the asymptotic value for T\ = 0.147 providing the good description 
of r^(Q 2 ) for Q 2 of the order of several GeV 2 , when the 3-loops radiative correction 
is included. This procedure may be considered as a sort of preliminary estimate, since 
the full 3-loops analysis is not available. 

To generalize our approach to the neutron case, we use the difference between 
proton and neutron instead of the neutron itself. Although it is possible, in principle, 
to construct a smooth interpolation for the functions g\ themselves |17j . it does not 
fit the suggested general argument on the linearity in /i^, since Jf _?1 (0) is proportional 
t° \^A,n ~ which is quadratic and, moreover, has an additional suppression due to 
the smallness of isoscalar anomalous magnetic moment. So we suggest the following 
parametrization for the isovector contribution of It{Q 2 ), namely I^~ n (Q 2 ), above the 
matching point, where again only the 1-loop term is presented explicitly 



jp-n 
1 T,pert 



(Q 2 ) = 0{Q 2 - Q\ 



2M 2 

2, ' rr p_ri (i 



Q 2 



a s {Q 2 
2tt 



-c, 



« O p » 



+ c n - 



« O n » 



Q 2 Q 2 

Here the transition value Ql may be determined by the continuity conditions in 
a similar way. We get the value Ql ~ 1.04GeV , which is not too far from that of 
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the proton case. Concerning IJfiQ 2 ), which is given by Eq. (J8|i. we have not neglected 
G%(Q 2 ) 3 and we have used its very recent determination |19] . The result of this 
calculation is very slightly modified compared to the case where one assumes G^(Q 2 ) = 
and all the subsequent results involving the neutron were obtained with a non-zero 
Gg(Q 2 ). The asymptotic value of Tl~ n = 0.21 is dictated by the Bjorken sum rule. 
We also took the same value << O n >>= O.llGeV 2 as in the proton case. The 
plot representing T^~ n (Q 2 ) is displayed on Fig. 2 and agrees well with the very recent 
experimental data |2Hj . 

This may be considered as an argument in favour of the general picture of power 
corrections obtained in QCD sum rules calculations [T3j, where the neutron correction 
is small. However the quantitative comparison with the calculations in the framework 
of the chiral soliton model ^H] would require a more detailed analysis. 

Now we have all the ingredients to turn to the behavior of neutron integral, which 
is simply obtained from the difference Y P X {Q 2 ) — T^ n (Q 2 ). It is shown on Fig. 3 and we 
notice that the strong oscillation around Q 2 = lGeV 2 , we had in the previous analysis, 
is no longer there. 




0.2 0.4 0.6 0.8 1 

Figure 2: Our prediction for Y\~ n (Q 2 ) versus Q 2 . The thick line is the new analysis 
to be compared with the thin line, which represents our previous approach without 
corrections. 



3 We thank G.Dodge for making this suggestion. 
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Figure 3: Our prediction for r™(Q 2 ) versus Q 2 . The thick line is the new analysis 
to be compared with the thin line, which represents our previous approach without 
corrections. 

We now have all the ingredients to investigate the deuteron integral. Note that in 
this case the generalization of the GDH sum rule may be naturally decomposed into 
two distinct regions. 

The first one is the region of large Q 2 , where nuclear binding effects can be disregard, 
so that the deuteron structure function is the simple additive sum of the proton and 
neutron ones. As a result, the Q 2 — > limit of this intermediate asymptotics is defined 
by the sum of the squares of proton and neutron anomalous magnetic momenta. 

When nuclear binding effects are taken into account one should get instead the 
square of the sum of these anomalous magnetic moments. As they are known to be 
rather close in magnitude and of different sign, the result should be therefore very 
small. 

The difference between these two regimes should be attributed [21] to the deuteron 
photodesintegration channel, which is supported by existing explicit calculations [22] 
in the case of real photons. For virtual photons, this allows to estimate the Q 2 value, 
where binding effects start to play a role, to be of the order of m 2 . The simplest way 
to implement this reasoning [21] is to use the following expression: 

i?(Q 2 ) = \[d{Q 2 dQ - Q 2 ) J Q *n2 + e (Q 2 - QDMiQ 2 ) + W) (19) 

Here we introduced the nuclear scale Qdo ~ rn^ and neglected the square of the 
deuteron anomalous magnetic moment. The prediction is shown in Fig. 4 and seems 
to be in good agreement with the preliminary JLAB data [23 . 

Let us finally discuss the role of the elastic contribution. It must be definitely 
included [21] if one uses the operator product expansion (OPE), which is the essen- 
tial tool in determining the power corrections. However, we use the OPE only above 
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Figure 4: Our prediction for Tf(Q 2 ) versus Q 2 . The thick line is the new analysis 
to be compared with the thin line, which represents our previous approach without 
corrections. 



matching point, where the elastic contribution is small. At the same time, below the 
matching point the object those Q 2 behaviour is studied may be considered as a sort 
of fracture function, where only a partial summation over final states, excluding the 
elastic one, is implied. It is this function which may reach the GDH value at Q 2 = 0. 
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